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Abstract. The MAGIC telescope observed the re- 
gion around the distant blazar 3C 66A for 54.2 hr 
in 2007 August-December. The observations re- 
sulted in the discovery of a 7-ray source cen- 
tered at celestial coordinates R.A. = 2 h 23 m 12 s 
and decl.= 43°0.'7 (MAGIC J0223+430), coincid- 
ing with the nearby radio galaxy 3C 66B. A 
possible association of the excess with the blazar 
3C 66A is discussed. The energy spectrum of 
MAGIC J0223+430 follows a power law with 
a normalization of (1.7 ± 0.3 st at±0.6 syst ) xlO -11 
TeV 1 cm 2 s 1 at 300 GeV and a photon index 
T = 3.10 ± 0.31 stat ±0.2 syst . 

Keywords: gamma rays: observations — BL Lac- 
ertae objects: individual (3C 66A) — galaxies: indi- 
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I. Introduction 

As of today, there are 26 known extragalactic very 
high energy (VHE, defined here as E > 100 GeV) 
7-ray sources. All of them are active galactic nuclei 
(AGNs) with relativistic jets. With the exception of the 
radio galaxy M 87 and CenA all detected sources are 
blazars, whose jets (characterized by a bulk Lorentz 
factor r ~ 20) point, within a small angle (6 ~ 1 /T), to 
the observer. The spectral energy distribution (logarithm 
of the observed energy density versus logarithm of the 
photon energy) of AGNs typically show a two-bump 
structure. Various models have been proposed for the 
origin of the high-frequency bump, the most popular 
being inverse Compton scattering of low-frequency pho- 
tons. There have been several suggestions for the origin 
of the low-frequency seed photons that are up-scattered 
to 7-ray energies: they may be produced within the 
jet by synchrotron radiation (synchrotron self-Compton 
or SSC mechanism [ 1 1) or come from outside the jet 
(external Compton or EC mechanism [2]). Relativistic 
effects boost the observed emission as the Doppler factor 
depends on the angle to the line of sight. In case the jet 
angle to the line of sight is large, models that depend less 
critically on beaming effects are needed [3|. The VHE 7- 
ray emission of AGNs might also be of hadronic origin 
through the emission from secondary electrons [4J. 

3C 66A and 3C 66B are two AGNs separated by just 
6' in the sky. 3C 66B is a large Fanaroff-Riley-I-type 
(FRI) radio galaxy, similar to M 87, with a redshift 
of 0.0215 0, whereas 3C 66 A is a distant blazar. 



As pointed aout in [6], the often referred redshift of 
0.444 for 3C 66A is uncertain. It is based on a single 
measurement of one emission line only Q, while in 
later observations no lines in the spectra of 3C 66A 
were reported O. Based on the marginally resolved host 
galaxy, a photometric redshift of ~ 0.321 was inferred 
|9l . In this paper we report the discovery of VHE 7- 
ray emission located 6/1 away from the blazar 3C 66A 
and coinciding with the radio galaxy 3C 66B in 2007. 
Detailed results and discussion can be found in |10|. 

II. Observations and Data Analysis 

3C 66A underwent an optical outburst in 2007 August, 
as monitored by the Tuorla blazar monitoring program. 
The outburst triggered VHE 7-ray observations of the 
source with the MAGIC telescope following the Target 
of Opportunity program, which resulted in discoveries 
of new VHE 7-ray sources in the past. 

MAGIC has a standard trigger threshold of 60 GeV, an 
angular resolution of ~ 0.°1 and an energy resolution 
above 150 GeV of - 25% (see EH for details). The 
MAGIC data analysis is described in detail in 0JJ, ifTTl . 

Data were taken in the false-source tracking (wobble) 
mode pointing alternatively to two different sky direc- 
tions, each at 24' distance from the 3C 66A catalog posi- 
tion. The zenith distance distribution of the data extends 
from 13° to 35°. Observations were made in 2007 Au- 
gust, September, and December and lasted 54.2 hr, out of 
which 45.3 hr passed the quality cuts based on the event 
rate after image cleaning. An additional cut removed the 
events with total charge less than 150 photoelectrons 
(corresponding to an energy of ~ 90 GeV) in order to 
assure a better background rejection. 

Just before the start of the observation campaign 
~ 5% of the mirrors on the telescope were replaced, 
worsening the optical point-spread function (PSF). As a 
consequence, a new calibration of the mirror alignment 
system became necessary, which took place within the 
observation campaign and improved the PSF again. The 
sigma of the Gaussian PSF (40% light containment) 
was measured to be 3/0 in 2007 August 12-14, 2/6 
in 2007 August 15-26 and 2/1 in 2007 September 
and December. To take this into account, data were 
analyzed separately for each period and the results were 
combined at the end of the analysis chain. However, the 
realignment resulted in a mispointing, which was taken 
care of by a new pointing model applied offline using 
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Fig. 1. Left plot: Significance map for 7-like events above 150GeV in the observed sky region. The green cross corresponds to the fitted 
maximum excess position of MAGIC J0223+403. The probability of the true source to be inside the green circles is 68.2%, 95.4% and 99.1% 
for the inner, middle and outer contour, respectively. The catalog positions of 3C 66A and 3C 66B are indicated by a white square and a black 
dot, respectively. Right plot: |Alpha| distribution after all cuts evaluated with respect to the position of MAGIC J0223+430. A 7-ray excess 
with a significance of 6.0cr is found, which corresponds to a post-trial significance of 5.4 a. 



starguider information |[T2l . Considering the additional 
uncertainty caused by the offline corrections, we esti- 
mate the systematic uncertainty of the pointing accuracy 
to be 2' [CD. 

III. Results 

Figure Q] (left plot) shows a significance map produced 
from the signal and background maps, both smoothed 
with a Gaussian of a — 6' (corresponding to the 7-PSF), 
for photon energies between 150 GeV and 1 TeV. Loose 
background rejection cuts are applied to keep a large 
number of gamma-like events. The center of gravity of 
the 7-ray emission is derived from Figure Q] by fitting a 
bell-shaped function of the form 



F(x, y) = A - exp 



{x-x) 2 + jy-y) 2 
2a 2 



(1) 



for which the distribution of the excess events is assumed 
to be rotationally symmetric, i.e., a x — a y = a. The 
fit yields reconstructed coordinates of the excess center 
of R.A. = 2 h 23 m 12 s and decl.= 43°0.'7. The detected 
excess, which we name MAGIC J0223+430, is 6/1 
away from the catalog position of 3C 66A, while the 
distance to 3C 66B is 1/1. We made a study to estimate 
statistical uncertainty of the reconstructed position. The 
probabilities are shown in Figure Q]by the green contours 
corresponding to 68.2%, 95.4%, and 99.7% for the inner, 
middle, and outer contour, respectively. Using this study 
we find that the measured excess coincides with the cata- 
log position of 3C 66B. The origin of the emission from 
3C 66A can be excluded with a probability of 95.6% 
(2.0 a). Linearly adding the systematic uncertainty of 
the pointing for this data set (2', see above), i.e., shifting 
the excess position by 2' toward the catalog position of 
3C 66A, the exclusion probability is 85.4% (1.5 a). 



To calculate the significance of the detection, an 
Alpha I distribution was produced (see Figure Q] left), 
where Alpha is the angle between the major axis of 
the shower image ellipse and the source position in 
the camera. A signal of 6.0 a significance (pre-trial, at 
the position of 3C 66B) and 5.4 it (post-trial, using 30 
independent trials) has been calculated. 

Figure |2] shows the light curve of MAGIC J0223+430 
together with the flux of 3C 66A in optical wavelengths. 
As we integrate over 7-ray events from a wide sky region 
(~ 0.07 deg 2 ), we cannot exclude that 3C 66A con- 
tributes to the measured signal. The integral flux above 
150GeV corresponds to (7.3 ± 1.5) x 10~ 12 cm~ 2 s" 1 
(2.2% of the Crab Nebula flux) and is the lowest ever 
detected by MAGIC. The 7-ray light curve is consistent 
with a constant flux within statistical errors. These 
errors, however, are large, and some variability of the 
signal cannot be excluded. 

For the energy spectrum of MAGIC J0223+430, loose 
cuts are made to keep the 7-ray acceptance high. The 
differential energy spectrum was unfolded and is shown 
in Figure [3] The spectrum can be well fitted by a power 
law which gives a differential flux (TeV -1 cm" 2 s" 1 ) 
of: 



dN 



dE dA dt 



= (1.7 ±0.3) x 10" n (£;/300GeV) 



-3.1±0.3 



(2) 

The quoted errors are statistical only. The systematic 
uncertainty is estimated to be 35% in the flux level and 
0.2 in the power law photon index 1111 . As we cannot 
exclude that 3C 66A contributes to the measured signal, 
the spectrum shown in Figure [3] represents a combined 
7-ray spectrum from the observed region. 
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Fig. 2. Upper panel: MAGIC integral flux above 150 GeV in bins of 3 days (except for periods where the sampling was coarser). The grey 
dashed line indicates the average 7-ray flux. Lower panel: Optical light curve of 3C 66A as measured by the KVA telescope. During the 
MAGIC observations 3C 66A was very bright at optical wavelengths. 



IV. Discussion and conclusions 

A VHE 7-ray source MAGIC J0223+430 was de- 
tected in 2007 August-December. Given the position 
of the excess measured by MAGIC above 150 GeV, the 
source of the 7-rays is most likely 3C 66B. The VHE 
7-ray flux was found to be on the level of 2.2% Crab 
Nebula flux and was constant during the observations. 
The differential spectrum of MAGIC J0223+430 has a 
photon spectral index ofr = 3.10±0.31 and extends 
up to - 2 TeV 

The VERITAS collaboration recently detected VHE 
gamma-ray emission from 3C 66A lfl41 . with most of the 
signal coming from a gamma-ray flare 2008 September- 
October. The measured integral flux above 200 GeV is 
6% of the Crab Nebula flux, and the energy spectrum 
is characterized by a soft power law with photon index 
T = 4.1 (see Figure [3). In view of this detection, we 
note that if 3C 66A was emitting 7-rays in 2007 August 
to December then its flux was at a significantly lower 
level than in 2008. 

The position of the source found by MAGIC is not 
spatially compatible with 3C 66 A at a level of 1.5a. 
In addition to this discrepancy in the source posi- 
tion, the measured spectrum extending up to ~ 2 TeV 
poses an additional drawback to the association of 
MAGIC J0223+430 with 3C 66A. Due to the energy- 
dependent absorption of VHE 7-rays with low-energy 
photons of the extragalactic background (EBL, lfT31l ). 
the VHE 7-ray flux of distant sources is significantly 



suppressed. We investigated the measured spectrum of 
MAGIC J0223+430 assuming a blazar-like emission 
mechanism, and following the prescription of 11161 . de- 
riving a redshift upper limit of the source to be z < 0.17 
(z < 0.24) under the assumption that the intrinsic energy 
spectrum cannot be harder than T = 1.5 (T = 0.666). 
This assumption of T > 1.5 is based on particle 
acceleration arguments ifTTll . and the fact that none of 
the sources in the EGRET energy band (not affected 
by the EBL) have shown a harder spectrum, a trend 
confirmed by Fermi-LAT in their first three months of 
observations. The latter assumption of T > 0.666 can be 
considered as an extreme case of the spectrum hardness, 
suggesting a monochromatic spectrum of electrons when 
interacting with a soft photon target field 1 1 8 ] @ The 
derived redshift upper limit is incompatible with the 
association of MAGIC J0223+430 with 3C 66A, unless 
the redshift of 3C 66A is much lower than the assumed 
value of z — 0.444, or exotic scenarios are responsible 
for its VHE emission ll20l . 

3C 66B is a FRI radio galaxy similar to M 87, which 
has been detected to emit VHE 7-rays fl2~Tl . Il22l . Since 
the distance of 3C 66B is 85.5 Mpc, its intrinsic VHE 
luminosity would be two to eight times higher than the 
one of M 87 (22.5 Mpc) given the reported variability 
of M 87 (221, El. 

As in the case of M 87, there would be several possi- 

2 See also 1191 for more detailed calculations. 
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Fig. 3. Differential energy spectrum of MAGIC J0223+430. No 
correction for the 7—7 attenuation due to the EBL has been made. The 
Crab Nebula spectrum ['11"] and the VERITAS spectrum on 3C 66A 
[14 ] is also shown as a reference (dashed grey line). 



bilities for the region responsible of the TeV radiation in 
3C 66B: the vicinity of the supermassive black hole E4l . 
the unresolved base of the jet (in analogy with blazar 
emission models; O) and the resolved jet. Unlike for 
M 87, we do not observe significant variability in the 
VHE 7-ray flux and therefore we have no constrains on 
the size of the emission region. However, as the angle to 
line of sight is even larger than in M 87 (M 87: 19°, ll25l : 
3C 66B: 45°, 11261 ) the resolved jet seems an unlikely 
site of the emission. On the other hand, the unresolved 
base of the jet seems a likely candidate for the emission 
site as it could point with a smaller angle to the line 
of sight. If the viewing angle was small, blazar-like 
emission mechanisms cannot be excluded. The orbital 
motion of 3C 66B shows evidence for a supermassive 
black hole binary (SMBHB) with a period of 1.05±0.03 
years [27 j. The SMBHB would likely cause the jet to 
be helical, and the pointing direction of the unresolved 
jet could differ significantly from the direction of the 
resolved jet. 

It has also been proposed that the spectrum obtained 
by MAGIC could be interpreted as combined emis- 
sion from 3C 66B, dominating above ~ 200 GeV, and 
3C 66 A l28l . The high-energy emission from the blazar 
would be strongly attenuated by the interaction with the 
EBL, and could only contribute to the measured excess 
at low energies. The emission from 3C 66B could be 
explained as radiation coming from the layer of the jet, 
in the framework of the structured jet model [29 1 already 
used to interpret the TeV emission of M 87 Q. 

Given the likely association of MAGIC J0223+430 
with 3C 66B, our detection would establish radio galax- 
ies as a new class of VHE 7-ray emitting sources. 
According to 1291 , there are eight FR I radio galaxies 
in the 3CR catalog that should have a higher 7-ray 
flux at lOOMeV than 3C 66B, but possibly many of 
these sources are rather weak in the VHE 7-ray band. 
Further observations of radio galaxies with the Fermi 



Gamma-ray Space Telescope as well as by ground-based 
telescopes are needed to further study the 7-ray emission 
properties of radio galaxies. 
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